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ABSTRACT 



The Boeing PAN AIR Technologies code (A502i) is investigated to explore its 
suitability for determination of separation forces on ordnance. To this end, A502i is first 
assessed by applying it to three problems for which other solutions and experimental data 
are available, i.e. steady flow past a rectangular, parabolic arc wing and a delta wing at 
both subsonic and supersonic conditions. Good agreement is found in all cases. A502i is 
then applied to the GBU-24’s being in two configurations for a subsonic case and a 
supersonic case. Good agreement is found with data obtained from wind tunnel 
experiments for low angles of attack. 
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I. INTRODUCTION 



In the past, ballistic trajectory determination for manual or computer predicted 
ordnance delivery from an aircraft was determined through measurement of separation 
forces on the piece of ordnance via wind tunnel or captive carry measurements. The advent 
of panel method codes using linearized potential theory, such as A502i, or its full potential 
version TranAir, allow for a cheaper and safer method of predicting separation forces. 
Furthermore, A502i allows for any arbitrary configuration to be modelled within the 
limitations of the number of panels and networks allowed and excluding transonic flow. 

The purpose of this work was to determine the separation forces on a GBU-24 
carried by an F-14 on stations 3 or 6 or both. It was also the purpose of this work to 
provide an analysis of the code itself to see if it is a viable tool for the study of flow 
characteristics over arbitrary wing configurations for use in the Naval Postgraduate 
School’s (NPS) Department of Aeronautics and Astronautics. The majority of the work 
was conducted on the NPS computer systems. The Department of Aeronautics and 
Astronautics Silicon Graphics Incorporated (SGI) workstations were utilized for most of 
the input files as well as the execution of the code. Due to the amount of disk space 
required, storage of the output files took place on the NPS Computer Center’s Y-MP EL98 
Cray computer. The bulk of the GBU-24 data was calculated using the SGI workstations 
at the Naval Air Warfare Center in Warminster. 

The scope of this analysis was to understand the capabilities of the A502i code. 
The approach was to validate A502i against existing data and linear theory. The code was 
run for three different geometries under assorted Mach and AOA conditions. Comparisons 
were made for each of the geometries. 
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II. OVERVIEW OF THE A502i CODE 



The A502i code is used to computationally analyze inviscid subsonic or supersonic 
flows about arbitrary configurations. The code differs from other panel methods in that it 
is a higher order panel method; that is, the singularity strengths are not constant on each 
panel. A502i solves the linearized potential flow boundary-value problem at subsonic and 
supersonic Mach numbers. 

The aerodynamic solution provides surface flow properties (flow directions, 
pressures, Mach number), configuration forces and moments, sectional forces and 
moments, and pressures. Additionally, A502i calculates flow properties in the flow-field 
points and flow-field streamlines. Results are limited to subsonic and supersonic cases 
(transonic cases excluded) with attached flow. Results are not usually applicable to cases 
where viscous effects and separation are dominant. 

A. THEORY 

The basic equations describing the flow of a viscous, compressible, heat- 
conducting fluid are the Navier-Stokes equations. These are: 

(a) The continuity equation, 



dp 

at 



+ V.(pV) = 



ap 

at 



+ L 



i-i 



a(pVj) 

a Xj 



(2.1) 



w here V = ( , - , ^ ^ - - ) is the gradient operator with respect to the location vector 

x = (x,,x,,x 3 ), and where the conventional index notation is used instead of x=(x,y,z). In 
addition, t is time, p(x,t) is the density and V(x,t) is the velocity vector, with components 
V = (V,,V 2 ,V 3 ). 
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(b) The conservation of momentum equation. 



IffpV^Z^fpV.Vp-^I^T. + pf, (j. 1.2.3) (2.2) 

where x ij is the deviatoric portion of the stress tensor which vanishes for a frictionless 

fluid, F(x,t) is an external body force per unit mass exerted on the fluid, and p(x,t) is the 
pressure. 



(c ) The conservation of energy equation 



d 1 ^ ^ i 

(pe + -j P I' V| 2 + P) + I Kpe + j p | V| : + p)V, ] 

dp „ d aT 

■dr +I ^ <T " V " tk aZ )+pl:fiVi 



(2.3) 



where e(x,t) is the internal energy of the fluid, k is the coefficient of heat conductivity for 
the fluid, and T(x,t) is the temperature. 

(d) The equation of state 



f(p,p,T)=0 (2.4) 

where the function f depends on the type of fluid. The derivations of these equations can 
be found in Refs. 1,2 and 3. 

The Navier-Stokes equations can be simplified by the neglect of viscosity, which is 
equivalent to setting the deviatoric stress tensor to zero. Combining the momentum and 
continuity equations yields 



dv : a p 

p_ dT = _ ax^ + pf j J = I ’ 2 ’ 3 

where the convective derivative operator is defined as 

— = — + y v -A_ 

dt at f 'ax, 



(2.5) 



Equation (2.5) is referred to as Euler’s equation. The continuity and energy equations can 
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Equation (2.5) is referred to as Euler's equation. The continuity and energy equations can 
be reduced to 



p-^ ( y|v| :) = -V. Vp +P V.? (2.6) 

and the rate of increase of heat per unit mass is given by 

1 dp d 1 

q ,_V.<kVT)--£ + p-(-) (2.7) 



Equations 2.5, 2.6 and 2.7 can be reduced to a single equation when four further 
assumptions are made. First, assume isentropic flow, thus 

q=0 (2.8) 

Second, assume irrotationality 



VxV = 0 (2.9) 

which allows for the introduction of the potential function [Refs 231 

VO> = V (2.10) 



Third, assume the existence of a freestream potential <1> b , whose gradient is the uniform 



velocity V m attained at points sufficiently distant from the disturbance being analyzed, and 
thus write 

<p = <J>- <|> (2.11) 

and 

V = (u,v,w) = Vd> = Vd> £ + V<j> = V,+ V(j) (2.12) 

The quantities <j> and v are called the perturbation potential and velocity [Ref 3[. Fourth, 
assume that 

|v| 2 <a, (2.13) 

everywhere, where a^ is the freestream speed of sound. 
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Based on these four assumptions, the unsteady potential equation is obtained 
| Refs 1,3]: 



( 1 - + <j> yy + <j> a - 2M 2 <j> Kt - M" tj> tt 

= M 2 [-^-(y- l)(2u + 2(j) t + 1 v |“ ) V 2 <j> 

+ (2u - u 2 )(j>„ + v 2 <j) yy + 2vw<j> yz + w 2 ^ 

+ 2(1 + u)(v<j) xy + w<]) s2 ) + 2(uu t + vv t + ww t )J 

(2.14) 

Assuming the flow conditions do not change with time yields the steady non-linear 
potential equation. 

(l -M 2 )4>„+<|) y y+<t> zz 
= M 2 [y(y- l)(2u+ |v| 2 ) V 2 <j) 

+ (2u+ u 2 )<[>„ + v : (j> yy + 2vw(j> y2 + w 2 ^ 

+ 2(1 + u)(v(j), y + w<t>„)] 



When M^O, equation (2.15) reduces to Laplace’s equation, 

V> = 0 



(2.15) 



(2.16) 



For the case of M, * 0, the following is supposed, 

M 2 |v|<1-M 2 (2.17) 

M 2 |v|« 1 (2.18) 

which are small perturbation assumptions [Refs. 1,2]. With these assumptions the steady 
non-linear potential equation reduces to the Prandtl-Glauert equation [Ref 1]: 



( 1 “ M 2 )<[>„ + <j)yy + <[> a = 0 



(2.19) 
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Through a coordinate transformation |Refs. 1 ,2,3 1, the Prandtl-Glauert equation can be re- 
written as: 

s<t>« + 4>„ + <l> a -0 (2.20) 

where when s=l, it is the subsonic case and Laplace’s equation applies and when s=-l, it 
is the supersonic case and the wave equation applies. Applying Green’s third identity [Ref. 
1 ] yields the following integral equation, 

<KP) = - — JJ[~nfi*vi-]dS (2.21) 

4k * k k 

where a represents the source strength and (j. represents the doublet strength. When 

supplemented with boundary conditions, it is equation (2.21) that A502i solves. 

A502i solves equation (2.21) through a discretization process. The general idea of 
the process falls into two parts. The first is developing finite dimensional approximate 
representation formulas for the singularity functions, which creates a system of equations 

with unknown coefficients, X ( . The second part involves solving the set of equations for 

all This allows for completely determining the source and doublet functions. Then, by 

virtue of equation (2.21), the potential function <}>(P) is determined for all points P, solving 
the problem. 

The features of A502i which distinguish it from predecessors are three-fold. The 
first is a feature known as “continuous geometry”, the second is linear source and quadratic 
doublet variation, the third is continuity of doublet strength. 

Most panel methods approximate the configuration geometry with panels whose 
planform is a quadrilateral. Thus, if the panels themselves are planar, only a small class of 
configurations (such as cylinders and flat wings) can be described without gaps being left 
between panels. These gaps are generally small, except for highly twisted surfaces. The 
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gaps cause little numerical error in subsonic flow, but in supersonic flow, the cumulative 
effect of the gaps is serious IRef. 1 ]. The problem is not associated with leakage of flow 
through the gaps, but with the doublet strength jumping abruptly from a non-zero value to 
zero at a panel edge which does not exactly meet the adjacent edge. In A502i, gaps are 
closed by means of panels which are comprised of several planar regions. 

The feature of linear source and quadratic doublet variation is what makes A502i a 
higher order panel method. The basis function corresponding to a source parameter is 
locally linear, while the basis function corresponding to a doublet parameter is locally 
quadratic. This is what allows for A502i to find supersonic solutions. Numerical solution 
of the wave equation is far more sensitive to the numerical idiosyncracies of a panel method 
than is the solution of Laplace’s equation for subsonic flow. Experimental evidence [Ref. 
1 ] indicates that exact surface analysis is not feasible in supersonic flow without doublet 
continuity, thus the potential for numerical error is greatly reduced by requiring the doublet 
singularity strength to be continuous across panels. 

B. GENERAL A502i USAGE 

The use of the A502i code consists of generating an input file, which can be 
arbitrarily named, and which contains the information defining the geometry of the 
configuration, flow field points of interest, the flow conditions and wakes. The process of 
building a geometry is difficult in that A502i is particular about its input format. Simple 
configurations, such as a rectangular, planform wing can be modelled manually, but more 
complex structures require a pre-processing program, such as MACGS, where a geometry 
can be graphically built. MACGS will output a data file in a format that, with minor 
modifications, via another pre-processing program that can move the data from three 
columns to six columns, will be readily usable by A502i. Currently, the school does not 
have a copy of MACGS, but it can be acquired through McDonnell-Douglas. To complete 
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this thesis, MACGS was used on the SGI workstations at NAWC Warminster. Wakes 
also must be constructed in the same manner as the structure to be analyzed. More detailed 
instructions on the specifics of wakes and surface geometries can be found in Ref. [4J. 
Appendix A is a portion of an output file, but lines 1 thru 1 120 are an exact duplicate of the 
input file. 

1. Running A502i with an Existing Executable 

Assuming an A502i executable file (e.g., A502) has already been placed in a user 
executable directory (e.g., /usr/local/bin), the only other necessary items needed to produce 
a set of A502i output files is the input file and a large amount of storage space. Anything 
modelled with more than one thousand panels total will use more than one hundred mb of 
disk space. If the maximum number of panels (20,000) is used, the disk space required 
will be on the order of 2 gb. 

To run A502i, enter after the UNIX prompt: 

A502 <input file > output file 

Prior to running the code, it is highly recommended that a Cray account be opened and 
linked to the department’s SGI workstations. This is done by assigning the same user i.d. 
number to the Cray account as is assigned to the account with the department. User i.d. 
numbers can be changed by the computer center at the user’s request. This is required due 
to the limited disk space available to individual accounts in the department. Once an 
account is opened, log on to a department workstation, change directories to an existing 
Cray directory, for example (after the UNIX prompt): 

cd /jedi/dl/maletour 

Transfer the input file to the Cray directory and execute the code. The screen will display 
what portion of the code it is performing and how long it took to perform each portion in 
CPU time. The code outputs numerous files in addition to the arbitrarily named output file. 
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The two output files of interest, in the vast majority of cases, are the arbitrarily named 
output file and the ft 13 file. In order to run another solution all output files must be deleted 
or renamed prior to re-executing the code. Relevant results should not be kept on the Cray 
account as files on disks dl, d2 and ul are considered temporary storage and subject to 
erasure after a period of time. 

2. Creating an A502i Input File 

The input file, which can be arbitrarily named, consists of two portions, the largest 
being the geometry data. Appendix A is a complete recreation of the input file for the 
GBU-24 with canards. The file begins with line 1, $TITLE, and ends with line 1120, 
$END. The line numbers are for reference only and are not part of the actual input file. 
The first portion consists of creating the initial conditions, i.e., the free-stream Mach 
number and angle of attack, the type of analysis to be performed, i.e., solution or 
datacheck, what types of output that are to be included in the arbitrarily named output file, 
and reference points to be used in calculating forces and moments. The geometry data 
consists of the points that bound each panel, that in turn belong to a specific group of 
panels that make up a network. The overall structure being modelled consists of a series of 
networks. A502i can run up to 150 networks and or 20,000 panels with a limit of 8,000 
panels per network. Referencing Appendix A, line 28 represents the first network of the 
model, a canard. Line 29 represents the number of networks that will be classified under 
this SPOINTS statement. Line 30 indicates what kind of surface the network will be, a 
three-dimensional surface with flow properties to be calculated, a wake and a base are 
several examples. Line 31 is the number of y points and the number of x points 
respectively that make up the grid of that network. Line 32 is where the panel points start. 
Reference 4 contains detailed instructions on the options and meaning of each of the non- 
geometry inputs, including some capabilities not shown in Appendix A. 
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Two types of solutions can be run, a datacheck and a full solution. Reference 4 
explains how to enter either one into the input file. The datacheck only analyzes the 
geometry. This can be accomplished in a matter of seconds for a simple geometry as it is 
only running the first several portions of the code. The full solution can take a couple of 
hours for a geometry of the size of 4,000 panels. The datacheck should be run once the 
geometry has been modelled. It will check for any panel edges that do not abut properly, 
and when column 4 of line 20 in Appendix A is a 1, the datacheck will list the unit normal 
vectors, which must be facing outward. The datacheck will also see if the wakes are 
attached properly. A502i is capable of giving warnings both on-screen and in the 
arbitrarily named output file when an edge or a wake is not modelled properly, but it only 
lists the unit normal vectors. The directions of the vectors must be manually checked by 
the user. The full solution performs the datacheck first, so the data is repeated in the 
arbitrarily named output file. Appendix B is a portion of the output file that contains the 
summary of facing surfaces. Each panel edge is looked at to see what it abuts against. 
Sections such as wingtips, leading edges of a flat plate or any surface that does not need a 
wake attached, but is unabutted to any other panel on that edge will draw probable error 
messages or warnings from the code. The user must ensure that the edge is not supposed 
to abut against anything or need a wake attached. If that is the case, the warnings may be 
ignored. Appendix D is the first page of the portion of the output file that lists the unit 
normal vectors. The three columns under zc are the x-y-z coordinates of the given panel’s 
center. The three columns under znc are the x-y-z coordinates of the unit normal vector. In 
most cases, when the y coordinates are of the same sign, then the unit normal vector is 
pointing outward. 



C. GEOMETRY MODELLING 



Five geometries needed to be modelled, each of increasing complexity. Modelling 
proved to be the most difficult task, in that A502i is a FORTRAN code and is very format 
sensitive, but the sheer number of points that need to be generated can take a lot of time and 
the order those points are listed in the input file is what determines whether or not the shape 
is correctly modelled. Of the five geometries modelled, none were done completely 
manually. A spreadsheet was used for generating the parabolic arc airfoil and the deltawing 
since those structures can be constructed out of one network, excluding wingtips and 
wakes, and the surface can be defined by a mathematical function. The bombs and the F- 
14 required the use of MACGS to be properly modelled. MACGS is indifferent as to the 
order that geometries are built, and often doesn’t require many coordinate inputs if building 
a model on top of an existing IGS file. The output file from MACGS is automatically 
formatted and the points placed in the appropriate order for A502i to understand. 
Although, the order may be reversed where the unit outward normal vector is concerned. 
MACGS has the ability to output files in several different panel method code input formats, 
including PMARC. Reference 4 gives detailed instructions on how to properly order 
points to build a group of networks that will model a geometry. A502i uses a right-handed 
coordinate system that is similar to an aircraft body axes. When put in terms of a wing, the 
x axis is positive from leading edge to trailing edge. The z axis is positive up and the y axis 
is positive out the right wing 

1. The Parabolic Arc Airfoil 

The parabolic arc airfoil is the simplest of all the geometries. The airfoil has a chord 
of five and a span of ten. The maximum thickness is .15. The model consists of 
approximately 600 panels, including the wake and wingtips. A spreadsheet was used to 
develop the geometry portion of the input file. Line 32 of Appendix A demonstrates the 



format that the spreadsheet used. Rows consist of two points, with coordinates xl , yl, zl, 
x2, y2, z2 using a format of 6F10.0. The chord was divided into 25 points (x coordinate) 
from trailing edge to leading edge and then another 25 points from leading edge to trailing 
edge (bottom half). The span was divided into 12 points (y coordinate) from left to right. 
Due to the symmetry of a rectangular planform, the y coordinate was constant along the 50 
x coordinates that constituted a chordwise cross-section. To attain a maximum thickness of 
. 15 the formula, 

z=.3 t) (2.22) 

c c" 

was utilized to generate the values of the z coordinates. The wingtips simply connect the x 
coordinate on the top side with it’s symmetrical counterpart on the bottom side. Due to a 
trailing edge composed of a straight line, the wake is modelled by a single panel that spans 
the trailing edge and has a length aft of 100. Figure 2. 1 shows the panel distribution across 
the top surface of the parabolic arc airfoil, where the thickness is represented by the color 
scheme. A panel and a point are numbered to show how they were entered into the input 
file. 

2. The Deltawing 

The deltawing represented a step up in complexity over the parabolic arc airfoil. 
The chordwise cross-section is parabolic, while the spanwise cross-section is linear. The 
procedure for building the geometry on a spreadsheet was the same as that for the parabolic 
arc airfoil, only the chord length is not constant along the span. For simplicity in design, 
the number of panels per column of panels is constant on the deltawing, as on the parabolic 
arc airfoil. This means an increasing panel density in the direction of the wing tip. The 
wake is modelled the same as the airfoil. The right wingtip ended in a point, so no extra 
panelling was needed to close any gaps. The symmetry of the deltawing allowed for 



13 




Figure 2.1 Parabolic Arc Panel Distribution 

further simplification and reduction of the code’s run time by only modelling from the 
centerline to the right tip. A502i allows the user to stipulate whether there is symmetry in 
the x-z plane and or the x-y plane (see line 5 and 6 of Appendix A). This means that the 
gap between the top and bottom panels at the center line does not need to be bridged as in 
the parabolic arc airfoil (symmetry could have also been used in the airfoil’s case). The 
chord of the deltawing has a length of 90 and the semi-span has a length of 15. The 
maximum thickness occurs midway along the centerline and is .05. The model consists of 
880 panels. Figure 2.2 shows the panel distribution along the top surface of the deltawing. 
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Thickness is represented by the color scheme. 




Figure 2.2 Panel Distribution of Deltawing 

3. GBU-24 

Wind tunnel experiments were run on the GBU-24 without canards attached, so it 
was deemed relevant to build a model with and without canards as a comparison of the 
code’s performance. The model of the GBU-24 was too complex to build with a 
spreadsheet, so MACGS was used. The bomb was modelled at NAWC Warminster by 
superimposing a group of networks on top of an IGS file being displayed by MACGS. 
The complete configuration consists of approximately 1300 panels. Figures 2.3 and 2.4 
are displays of the GBU-24 with canards, with Figure 2.4 including the wakes. Figure 2.5 
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is included to show how the GBU-24 model was assembled. Each different color 
represents a network. 

Several features of the geometry are relevant to point out. Two of them are 
modifications made to the geometry that differ from the actual dimensions of the bomb. 
Dr. Alex Cenko of NAWC Warminster has extensive experience with modelling stores in 
A502i. The modifications were made on his knowledge of how to get the most accurate 
results from the code when modelling stores. The first is to model the fins and canards as 
flat plates, i.e., no thickness, which A502i allows you to do through a single numerical 
change in the input code for each network that represents a flat plate (see line 30 of 
Appendix A). The fins and canards are extremely thin when compared to the rest of the 
bomb, and to add a third dimension to the geometry complicates the construction of the fin 
or canard for several reasons. The leading and trailing edges must be sharp and the surface 
the fin or canard attaches to would have to be modified to abut properly with two edges 
instead of one. Experience has shown that the simpler version yields accurate predictions. 
A502i is an inviscid code, so it cannot take into account separation effects on its own. The 
GBU-24 does not have a flat base. In reality, it is more bullet nosed in shape. However, 
at the speeds with which the bomb is being analyzed, separation does occur near the trailing 
edge of the bomb. Experience has shown that truncating the end into a flat base and 
designating it a separated flow region through an appropriate input (see line 702, column 1, 
Appendix A) yields better results than attempting to model the bomb to exact physical 
dimensions. The last feature to point out are the wakes, as seen in Figure 2.4. A502i has a 
limitation in that the wakes must be modelled by the user, and they have the same abutment 
requirements as physical surfaces. Regardless of angle of attack, the wakes remain 
stationary with respect to the body to which they are attached. At higher angles of attack, 
the wakes are no longer close to paralleling the free-stream velocity. Remodelling the 
wakes is nearly an impossible task. The fin and base wakes would not be too difficult 
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Figure 2.3 GBU-24 Geometry (Wakes not Shown) 
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Figure 2.4 GBU-24 Geometry (Wakes Shown) 
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Figure 2.5 GBU-24 Geometry (Network Modelling) 
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because they do not abut against any physical surfaces except for the surfaces they are 
trailing from. The canard wakes must abut properly against the bomb's body all the way to 
the base. Modification of those wakes would entail modification of the entire body, or in a 
possible simplification, letting the wakes remain attached along the body until the base of 
the bomb and then shifting them relative to the free-stream. 

4. The F-14 

The F-14 geometry was modelled in the same fashion as the GBU-24. The 
geometry consists of approximately 1500 panels. While that may seem fairly coarse for 
such a complex structure, experience shows that it is all that is required to get accurate 
predictions. The primary area of interest is the underside of the fuselage forward and 
between the two nacelles. Higher panel density on the top half is not required. Figures 2.6 
and 2.7 display the F-14 geometry without and with wakes shown. Several omissions are 
made to the model as having a trivial effect on the analysis or no effect at all. Phoenix rails 
and bomb racks are not modelled along with the chin pod because they are deemed 
insignificant to achieve reasonably accurate predictions over small angles of attack. 
External tanks were not considered, but could be modelled much in the same way as the 
bomb and inserted into the input file to see what effects the drop tanks have on separation 
forces. The vertical tails and horizontal stabilizers were deemed irrelevant to the prediction 
of the separation forces and were left out. This reduces the number of panels and 
networks, which also reduces the amount of time it takes to run a solution. 

5. Combination Geometries 

The F-14 and GBU-24 were modelled separately, but were combined together as 
shown in Figure 2.8. The first step to accomplish this was using the FORTRAN code 
NAVSEP which, among many of its functions will translate coordinates to relocate items in 
the flow-field. Once accomplished, the GBU-24 file was pasted into the F- 14 input file. 



20 



Figure 2.6 F- 14 Geometry ( Wakes not Shown) 
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Figure 2.7 F-14 Geometry (Wakes Shown) 
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Figure 2.8 F- 14 and GBU-24 on Station 3 Combined Geometries 
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D. GRAPHICS VISUALIZATION 



One of the pre/post-processing codes that came with A502i is called RAID. 
Currently an executable exists in the department’s computer system that can be accessed 
through typing after the UNIX prompt, 

raid 

Raid is a basic graphics program that can read A502i and TRANAIR input files and display 
geometries and flow properties from solutions. After accessing raid, it will ask what type 
of input file it is being asked to visualize. It can handle five other modes of input besides 
A502i/TRANAIR [Ref. 10]. It will then ask for the name of the input file. A prompt will 
follow asking about object definition matrices which needs to be answered by, 

EACH 

The next prompt will ask if the panels are going to be shaded by a Cp value. Cp is used 
generically in that Cp can be displayed or Mach or any of the relevant 49 surface flow 
properties [Ref. 4]. This is only used after a solution has been run, data has been extracted 
and a colorscale for contour plots has been determined. If only the geometry is to be 
displayed, then hit carriage return for all the next questions until a pink window appears 
with a menu in the lower left comer of the window. An anomaly of the program is that if 
you want to display the wakes, then wake display must be deselected. To view the 
geometry select view on the menu. All selections with the mouse in RAID are made with 
the center mouse button. A new window appears with a menu bar at the bottom and left 
and the geometry in the center. From there, rotation, translation, scaling, axes, reflections 
and other manipulation of the geometry is possible. Figures 2.3 through 2.8 are examples 
of geometries displayed on RAID. 

When presentation of solutions (i.e., Cp or Mach contours) are desired, the use of 
another post-processing program is required to generate a colormap file. The program is 
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called crebar. An executable currently exists on the department's system. Type in, 

crebar 

after the UNIX prompt. The program will ask straightforward questions. Number values 
associated with colors available can be found in Ref. 10. The color file can be saved under 
any name, but must lie in the same directory as the input and solution files. The first line of 
the color file will list four numbers. The first number is the number of colors assigned to 
the colormap (248 maximum). Occasionally, the color bar displayed in RAID when using 
a colormap will disappear when certain menu items are selected. To prevent this, change 
the last three numbers to read 6, 1,-1. Plotting outputs from RAID requires saving the file 
in a format, such as RGB, that a printer will recognize. It is possible to change the text 
color and background color, the default is black, to avoid excessive use of black ink in hard 
copies. 
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III. DATA EXTRACTION 



When a solution has run to completion, there are two files of interest, the arbitrarily 
named output file and the ft 13 file. The arbitrarily named output file contains results for 
everything that A502i solves for. The ft 13 file contains only the 49 surface flow properties 
on each panel. Appendix D is the solution portion of an arbitrarily named output file for the 
first network. 

For purposes of displaying flow properties on RAID, it is necessary to utilize the 
ft 13 file. A post-processing code called RA1DCONV is used to extract the specific 
information. To access RAIDCONV, type 

raidconv 

after the UNIX prompt. An executable currently exists in the department’s system. The 
ft 13 file must be in the same directory as RAIDCONV is accessed in. RAIDCONV will 
prompt the user for which kind of panel method is being used (A502i is one of three 
choices). The next prompt will ask for the name of the ftl3 file. The last prompt will ask 
for the flow property that is to be extracted. A file called ftl3.cp will be created. It can be 
renamed for purposes of multiple flow properties extraction. Abbreviations for the 49 flow 
properties can be found at the bottom of page 1 of Appendix D. The two primary flow 
properties are 

LMACHU for local Mach number 
CP2ND for second order pressure coefficient 
CP2ND is the default setting for RAIDCONV. Once the ftl3.cp file is created, RAID can 
be used as previously discussed to display the flow properties. An anomaly of 
RAIDCONV is that it does not recognize kt=20 type wakes, used where wakes from a 
wing abut against a body [Ref. 4|. To assist in extracting all the data, the kt=20 wakes 
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should be placed at the end of the input file. In general, it is good practice to place all 
wakes at the end of the input file when using A502i. 

The arbitrarily named output file duplicates the data found in the ft 13 file and 
includes moments and forces. A502i will sum up the moments and forces on each network 
and for all networks so far [App. D]. The moments are computed based on the coordinates 
entered into the input file [Ref. 4 and App. A|. 
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IV. RESULTS OF A502i COMPUTATIONS 



A. PARABOLIC ARC AIRFOIL DISCUSSION 



This simple geometry was analyzed primarily to evaluate A502i's capabilities by a 
comparison to known linear theory. To this end, the geometry discussed in section II-C 
and shown in Figure 2. 1 was run by A502i at a Mach of 0.3 and a Mach of 1 .5 at an angle 
of attack of zero. Two of the 49 flow properties that A502i computes [Ref. 4] for each 
panel are linear Cp and second order Cp, given by 

CPLIN = -2u c (4.1) 

CP2ND = -2uc-[( l-M*2) Uc 2+v c 2+w c 2] (4.2) 

Where u c , v c and w c are the compressible components of the perturbation velocity. 
Figure 4.1 plots the linear theory, A502i linear and second order results for the subsonic 
case, while Figure 4.2 represents the supersonic solution. 

Linear theory for parabolic arc airfoils is outlined in Refs. 2 and 3. The equation 
representing the subsonic case is given by: 



Cp(x) = 



-8 *x„ 



jt * chord * ^ 1 - M^ 



= *(1 -(.5-x)*ln 



1 — x 



) where 0< x< 1 



(4.3) 



Equation 4.3 includes a Prandtl-Glauert compressibility correction. The equation 
representing the supersonic case is given by. 



Cp(x) = 



20 



Vm;-i 



where 0 < x < chord 



(4.4) 



and 6 = x * ( 



l 



2x 



ma * chord chord 



-) 
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Figure 4.1 Cp Comparison of a Parabolic Arc Airfoil at Mach = 0.3. 




Figure 4.2 Cp Comparison of a Parabolic Arc Airfoil at Mach = 1.5. 
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Figure 4.1 shows very good agreement with the linear theory curve. There is a 
small but noticeable difference between A502i's linear results and second order results, 
with the second order results being more accurate, as expected. The maximum difference 
between the linear theory curve and the second order A502i curve amounts to a value of 
2.5% right at mid-chord. The gap between the two curves from .25 chord to .75 chord 
were the result of thin panel density in that region. 

Figure 4.2 shows excellent agreement with the linear theory curve. The A502i 
values of Cp for the linear and second order analysis are virtually identical. There are small 
deviations from linear theory near the leading and trailing edges, but this is expected due to 
numerical error associated with the discontinuity A502i would encounter right on the 
leading or trailing edges. 

B. DELTAWING DISCUSSION 

Reference 9 provides data on Mach distribution, using approximated linear theory, 
over a deltawing of the configuration discussed in section II-D. This simple geometry 
provided another test of A502i’s capabilities. Figures 4.3 and 4.5 show the A502i results 
for the subsonic and supersonic case, while Figures 4.4 and 4.6 reflect the results from 
Ref. 9. For both cases, good agreement is found with the linear theory, with A502i's 
subsonic analysis being physically more accurate than the approximate linear theory, while 
A502i’s supersonic analysis is not as physically accurate. 

A comparison of Figures 4.3 and 4.4 reveals several points of interest. The Mach 
contour representing the free-stream value is given by the dashed line. All lines outside the 
dashed line represent areas where the Mach value is less than free-stream, and inside the 
dashed line is where the Mach value is more than free-stream. The location of where the 
free-stream Mach contour, in Figure 4.3, intersects the centerline agrees very well with 
Figure 4.4. However, Figure 4.4 does not have the contour extending all the way to the 
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tip. This is a physical limitation of the approximate theory used in Figure 4.4 and A502i is 
giving a more realistic solution. Figure 4.4 suggests that the peak Mach value occurs at 
approximately two-thirds chord along the centerline. The A502i results show the peak 
Mach contour occurring out midway along the semi-span. Those Mach values are less than 
1% larger than the yellow Mach contour surrounding it, and can be attributed to how the 
panel density increases with movement towards the wingtip. A502i performed very well 
for this subsonic case. 

A comparison of Figures 4.5 and 4.6 shows that A502i did not perform as well as 
in the subsonic case. Again, the contour representing the free-stream value of 
Mach =1.414 is given by the dashed line. All lines forward are below free-stream and all 
lines aft are above free-stream. Figure 4.6 shows the intersection of the free-stream Mach 
contour on the centerline occurring at approximately 39% chord, which is in excellent 
agreement with A502i’s result. Figure 4.6 shows the peak Mach value occurring at the 
trailing edge on the centerline. This makes more physical sense than the results 
that A502i yielded. The maximum thickness of the deltawing occurs along the centerline, 
allowing for greater expansion. The discrepancy may be attributable to panel density and 
accumulation of numerical error. A close study of the A502i results reveals some 
discontinuities along the column of panels out at the wing-tip which would have adversely 
affected the solution and caused errors to propagate along the semi-span. 
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b> h> Chord 



Lines of Constant Mach (Minf=.01) 
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Chord 



Lines of Constant Mach (Minf=1 .414) 




Figure 4.5 A502i Mach Contour Plot (Moo=1.414) 





Figure 4.6 Approximate Linear Theory [Ref. 9] Mco=1.414 
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C. GBU-24 FREE-STREAM (NO CANARDS) DISCUSSION 



A free-stream measurement of separation forces on the GBU-24, without canards, 
was conducted in a wind tunnel for various Mach numbers from .8 to 1.2 |Ref. 7], Since 
A502i uses linear potential theory, the model of the GBU-24, without canards, was 
evaluated at both Mach .8 and 1.2, avoiding the transonic regime, to ascertain the accuracy 
of the code with the given geometry. Normal forces and pitching moments for both cases 
are plotted and compared to the wind tunnel data. 

1. Subsonic Case (Moo=0.8) 

The GBU-24 model, without the canards, was run for angles of attack varying 
from -10 to +10 degrees in two degree increments. Values much higher than that ran into 
wake modelling problems as the wake’s angle relative to the free-stream was getting large 
enough that results would become questionable, and remodelling the wake was too difficult 
for such a complex geometry. The results of the A502i analysis are displayed in Figures 
4.7 and 4.8. For angles of attack between -4 and +4 degrees, A502i does a good job of 
predicting the separation forces. The pitching moment, which happens to be unstable 
without the canards, is approximately linear over the -4 to +4 degree range and is the 
limiting factor to the models accuracy. The normal force is approximately linear over a 
wider range, and A502i does a good job of predicting the normal forces from -6 to +6 
degrees. 

2. Supersonic Case (M 00 — 1 • 2) 

The results for the subsonic case demonstrated that the effective range of angle of 
attack that A502i needed to explore was from -6 to +6 degrees. Figures 4.9 and 4.10 plot 
the comparison of wind tunnel data versus A502i results for pitching moment and normal 
force. The results for the supersonic case are slightly better than that of the subsonic case. 
The actual pitching moment of the GBU-24 is approximately linear over a wider angle of 
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GBU-24 Free-Sream (M=0.8) 
Canards Off 
Pitching Moment 




Figure 4.7 Comparison of Pitching Moments 



GBU-24 Free-Stream (M=0.8) 



Canards Off 
Normal Force 




Figure 4.8 Comparison of Normal Forces 
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attack region, but fluctuations in the data at -6 and +6 degrees means that the model is still 
only viable from -4 to +4 degrees. The normal force line is nearly linear from -10 to +10 
degrees and extrapolating the A502i results out to 10 degrees would still yield good 
predictions. 



GBU-24 Free-Sream (M= 1.2) 
Canards Off 
Pitching Moment 




Figure 4.9 Comparison of Pitching Moments 
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GBU-24 Free-Stream (M=1.2) 



Canards Off 
Normal Force 




Figure 4.10 Comparison of Normal Forces 

D. GBU-24 FREE-STREAM (WITH CANARDS) DISCUSSION 

As in the case with no canards, free-stream measurements of the separation forces 
on GBU-24 were taken from Mach .8 to 1.2 in a wind tunnel [Ref. 7|. Again, due to the 
limitations of linear theory inherent in the code, an analysis was done for Mach numbers of 
.8 and 1.2 to minimize transonic effects. Even with the more complex geometry, A502i 
does an accurate job of predicting the separation forces over the range of angles of attack 
that are approximately linear. 

1. Subsonic Case (Moo=0.8) 

The GBU-24 model, with canards, was run in two degree increments of angle of 
attack from -10 to +10. The wake modelling limitation, as well as reviewing the data from 
the wind tunnel measurements [Ref. 7] showed the non-linearity of the separation forces at 
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the higher values of angle of attack, precluded any attempts to predict forces beyond the 
aforementioned angle of attack interval. The results of the A502i analysis are displayed in 
Figures 4.1 1 and 4. 12. The addition of the canards makes the pitching moment stable, but 
linear over a smaller region than without the canards. A502i gave accurate results from -2 
to +2 degrees angle of attack when predicting pitching moment. The prediction of normal 
forces fared better, showing accurate results from -3 to +3 degrees angle of attack. 



GBU-24 Free-Stream (M=0.8) 
Canards On 




Figure 4.11 Comparison of Pitching Moments 
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GBU-24 Free-Stream (M=0.8) 



Canards On 
Normal Force 




Figure 4.12 Comparison of Normal Forces 



2. Supersonic Case (M<»=1.2) 

As in the case with no canards, the region of accuracy, with the model used, was 
assumed to be less than + or - 10 degrees angle of attack. Cases were run from -6 to +6 
degrees angle of attack in two degree increments. A comparison of A502i results with 
wind tunnel data is shown in Figures 4. 13 and 4.14. For both the pitching moment and the 
normal force, A502i does a much better job of prediction then when subsonic. The wind 
tunnel data is nearly linear in both pitch moment and normal force from -8 to +8 degrees 
angle of attack. Extrapolating the A502i data out to + or - 8 degrees angle of attack, shows 
excellent agreement with the wind tunnel data. 
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GBU-24 Free-Stream (M=1.2) 



Canards On 
Normal Force 




Figure 4.14 Comparison of Normal Forces 
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E. F-14 DISCUSSION 



The GBU-24, without canards was located at station 3 of the F-14 via the NAVSEP 
code. In the process of running a solution, the combined geometries were found to have a 
total of 150 networks (A502i's maximum). The combined geometry had about 4,000 
panels, far short of the 20,000 panel maximum, so there was room for more detailed 
modelling, but there was not a chance to insert the GBU-24, with canards, into the F-14’s 
flow field and analyze the forces on the bomb. The geometry with the canardless bomb ran 
to a solution that appeared to be valid, so there is a high degree of confidence that if the 
number of networks could be reduced to allow the GBU-24, with canards, to be inserted 
into the F-14 flow field, the code would yield accurate predictions at small angles of attack 
on the forward stations. To reduce the total number of networks by combining existing 
networks would have required a large time investment and the use of MACGS, which the 
department currently does not possess, the two reasons why it was not done. Figure 4. 15 
shows a Mach distribution of the solution of the canardless bomb and F-14 at Mach = 0.8 
and 0 degrees angle of attack 

F. POST PROCESSING DISCUSSION 

The Mach values for the subsonic case of the GBU-24, with canards, at 4 degrees 
angle of attack, were extracted from the ft 13 file. These values, used in conjunction with 
RAID are shown in Figures 4.16 and 4.17. The color distribution over the nose in Figure 
4.16 indicates that the bomb is at an angle of attack, and scanning the rest of the model 
showed no discontinuous solutions, which is generally represented in A502i by a Mach 
value of 0 or 1,000. The visual representation is a quick way of telling if A502i ran an 
accurate solution. The only other way is to individually check the Mach or Cp values of 
each panel in the ft 13 file or the arbitrarily named output file. The other point of interest in 
Figure 4. 16 is the lack of panel density along the mid-section of the bomb. The goal, in the 
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case of stores separation prediction, is to have as simple a model as possible that still gives 
accurate predictions. The fewer the number of panels, the shorter the run time. The fact 
that A502i is a higher order panel method allows the luxury of using fewer panels. Figure 
4.17 highlights the approach used to take into account separation effects as discussed on 
page 16. 
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Figure 4,15 Mach Distribution over GBU-24 and F-14 (Moo=0.8, aipha=0°) 
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Figure 4.16 Mach Distribution over GBU-24 (Moo=0.8, alpha=4°) 
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Figure 4.17 Mach Distribution over GBU-24 (Moo=0.8, alpha=4°) 
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V. SUMMARY AND CONCLUSIONS 



The main goal of this analysis is to determine the accuracy of A502i on both simple 
geometries and complex geometries. To accomplish this, A502i is compared with results 
available from linear theory and wind tunnel experiments. This allows conclusions to be 
drawn on the capabilities as well as the limitations of A502i. 

In general, A502i can accurately predict flow properties, forces and moments on 
simple and complex geometries at low angles of attack. The predictions are valid over a 
wide Mach range, from 0 up to and including 0.8 and from 1.2 and above. The supersonic 
solutions are available due to A502i’s higher order capabilities. 

The limitations of A502i are consistent with most panel methods. A502i cannot 
predict flow dominated by viscous, separated or transonic effects. It cannot predict flow 
with different total pressures, such as flow properties inside a jet plume or a propeller 
slipstream swirl. The biggest shortcoming of A502i is its inability to handle unsteady cases 
and automatically determine wake shapes. 

Experience or knowledge of the flow properties around the geometry being tested is 
important in building an accurate model. An accurate model may not be physically 
accurate. Flight test results revealed a yawing moment on the GBU-24 that was not 
discovered in wind tunnel experiments when the bomb was carried on an aft station. The 
yawing moment may be caused by the fact that the canards are not fixed, but spring- 
damped. A502i predictions would not be accurate without inputting a moment to simulate 
the deflection of the canards, since the canards are fixed by the geometry. A502i, while a 
powerful tool in terms of cost savings and time, cannot completely substitute for wind 
tunnel experiments and flight tests, as constructing a complex geometry to exact physical 
specifications will probably not yield accurate predictions. 
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Printed by maletour from osprey 



APPENDIX A. GBU-24 OUPUT FILE (INPUT FILE PORTION) 
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Printed by maletour from osprey 



Printed by maletour from osprey 
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183 23.730 - 3.138 2.324 

184 45.183 - 4.642 3.583 42.732 - 4.642 3.582 

185 40.281 - 4.642 3.582 37.830 - 4.642 3.582 

186 35.380 - 4.642 3.581 32.929 - 4.642 3.581 

187 30 478 - 4.642 3.581 28.027 - 4.642 3.580 
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Printed by maletour from osprey 
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Printed by maletour from osprey 



Printed by maletour from osprey 
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Printed by maletour from osprey 



Printed by maletour from osprey 
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Printed by maletour from osprey 



Printed by maletour from osprey 
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